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We investigate the chiral phase transition at nonzero temperature T and 
baryon-chemical potential /x^ within the framework of the linear sigma model 
pH I and the Nambu-Jona-Lasinio model. For small bare quark masses we find in 



both models a smooth crossover transition for nonzero T and fj,B = and a 



5h ' first order transition for T = and nonzero jj-b- We calculate explicitly the 

first order phase transition line and spinodal lines in the (T, plane. As 
expected they all end in a critical point. We find that, in the linear sigma 
model, the sigma mass goes to zero at the critical point. This is in contrast 
to the NJL model, where the sigma mass, as defined in the random phase 
approximation, does not vanish. We also compute the adiabatic lines in the 
(r, /is) plane. Within the models studied here, the critical point does not 
serve as a "focusing" point in the adiabatic expansion. 
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I. INTRODUCTION 



Chiral symmetry is spontaneously broken in the QCD vacuum. Lattice QCD simulations 
at nonzero temperature T and zero baryon-chemical potential indicate that chiral sym- 
metry is restored above a temperature T ~ 150 MeV [0. Such temperatures are believed to 
be created in nuclear collisions at ultra-relativistic energies. Consequently, a phase where 
chiral symmetry is transiently restored may be formed in these collisions. The subsequent 
expansion cools the system and takes it to the final hadronic state, where chiral symmetry 
is again spontaneously broken. 

It is important to determine the order of the chiral transition, as this infiuences the 
dynamical evolution of the system. For instance, it has been shown that a first order 
transition may lead to a defiagration wave and to a "stall" in the expansion of the system 
0). It has also been shown that a first order transition in rapidly expanding matter may 
manifest itself by strong non-statistical fiuctuations due to droplet formation . In the case 
of strong supercooling it may lead to large fiuctuations due to spinodal decomposition . 
In a second order phase transition one may expect the appearance of critical fiuctuations 
due to a large correlation length [Q. Experimentally, large-acceptance detectors are now 
able to measure average as well as event-by-event observables, which in principle allow to 
distinguish between scenarios with a first order, a second order, or merely a crossover type 
of a phase transition. 

Theoretically, the QCD phase diagram in the (T, /i^) plane has recently received much 
attention (see P-^). QCD with Nf = 2 fiavors of massless quarks has a global SU{2)l x 
SU{2)ji symmetry. This symmetry is spontaneously broken in the QCD vacuum, such that 
the order parameter (p^^ ~ {q^l^r) acquires a non-vanishing expectation value, where is 
the quark field are the fiavor indices). At zero baryon-chemical potential, the effective 
theory for this order parameter is the same as the 0(4) model which has a second order 
phase transition. Therefore, by universality arguments |10|, the chiral transition in Nj = 2 
QCD is likely to be of second order at /x^ = 0. Nonzero quark masses introduce a term 
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in the QCD Lagrangian which exphcitly breaks chiral symmetry. Then, the second order 
transition becomes crossover. 

At nonzero baryon-chemical potential, it is more difficult to infer the order of the chiral 
transition from universality arguments [|lT|. One commonly resorts to phenomenological 
models to describe the chiral transition in this case. Depending on the parameters of these 
models, they predict a ffist order, a second order, or a crossover transition. However, if there 
is a second order phase transition for /x^ = and nonzero T and a ffist order transition 
for small T and nonzero fiB, then there exists a tricritical point in the {T,^b) plane where 
the line of ffist order phase transitions meets the line of second order phase transitions. For 
nonzero quark masses, this tricritical point becomes a critical point. 

It has recently been proposed that this point could lead to interesting signatures 
in heavy-ion collisions at intermediate energies, if the evolution went through or close to 
this critical point. At this point, susceptibilities (e.g. the heat capacity) diverge, and the 
order parameter field becomes massless and consequently fluctuates strongly, which could 
be detected in event-by-event observables. 

In this paper we investigate the thermodynamics of two popular models of chiral dy- 
namics, the linear sigma model coupled to quarks , and the Nambu-Jona-Lasinio (N JL) 



model . Both models are tuned to reproduce correctly properties of the physical vacuum. 
Our goal is to study the chiral transition and to verify the existence of the critical point at 
nonzero chemical potential and temperature. We also study the behavior of isentropes in 
the vicinity of the phase transition line in the (T, /i^) plane. These results can then be used 
in dynamical simulations to confront the predictions of with experimental data. 

The structure of the paper is as follows. In Section II we study the thermodynamics of 
the linear sigma model coupled to quarks. This part of the paper is an extension of our 
previous study in ref. In Section III we do the same for the NJL model. Section IV 
presents numerical results. We conclude in Section V with a summary of our results. Our 
units are h = c = ks = I, the metric tensor is g'^^ = diag(+, — , — , — ). 
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II. THERMODYNAMICS OF THE LINEAR SIGMA MODEL 



The Lagrangian of the hnear sigma model with quark degrees of freedom reads 

C^q [i^^d^ - g{a + i^,T + \ {d.ad^'a + 9^7? • d^T^) - U{a, tt) , (1) 
where the potential is 

U{a, 7?) = ^ {a^ + 7r^-v^y-Ha. (2) 

Here q is the light quark field q = {u,d). The scalar field a and the pion field tt = (tti, 7r2, tts) 
together form a chiral field $ = (c, tt). This Lagrangian is invariant under chiral SU{2)l x 
SU{2)fi transformations if the explicit symmetry breaking term Ha is zero. The parameters 
of the Lagrangian are usually chosen such that the chiral symmetry is spontaneously broken 
in the vacuum and the expectation values of the meson fields are (a) — and (tt) = 0, where 
— 93 MeV is the pion decay constant. The constant H is fixed by the PC AC relation which 
gives H = f^Km^i where m^r = 138 MeV is the pion mass. Then one finds = — m^/)?. 
The coupling constant is determined by the sigma mass, = 2A^/^ + m^, which we set 
to 600 MeV, yielding ^ 20. The couphng constant g is usually fixed by the requirement 
that the constituent quark mass in vacuum, Myac = 5 A, is about 1/3 of the nucleon mass, 
which gives gr ~ 3.3. 

Let us consider a system of quarks and antiquarks in thermodynamical equilibrium at 
temperature T and quark chemical potential fi = iib/'^- The grand partition function reads: 

Z = Tt exp[- (n- ^Af^ /t] = J VqVqVaVjr exp j (C + q-^"^ q) . (3) 

Here = i j^^^ dt Jy d^x, where V is the volume of the system. We adopt the mean-field 
approximation, replacing a and tt in the exponent by their expectation values. Then, up to 
an overall normalization factor: 

Z = exp (--^) / 'DqVq exp |^ q [i'j^df, - g{a + i'jsT ■7r)]q + //g7°g| 

- exp (-^) detp { [p^Y + /X7° - g{a + i^^r ■ tt)] /t} . (4) 



All thermodynamical quantities can be obtained from the grand canonical potential 

rj-i 1^ 

Q{T, fi) = ^ = U{a, vr) + Q,^, (5) 



where the quark and antiquark contribution reads: 



Q,,-(T, /x) = -u, I |i? + T In [l + exp (^^)] + T In [l + exp [-AjT^)] } • (6) 

Here, z/g = 2NcNf = 12 is the number of internal degrees of freedom of the quarks, Nc = 3, 
and E = \J'p^ + M'^ is the valence quark and antiquark energy. The constituent quark 
(antiquark) mass, M, is defined to be: 

= g\a'' + TT^) . (7) 

The divergent first term in eq. (^) comes from the negative energy states of the Dirac 
sea. As follows from the standard renormalization procedure it can be partly absorbed in 
the coupling constant and the constant . However, a logarithmic correction from the 
renormalization scale remains and is neglected in the following. Similar logarithmic terms 
are explicitly included in calculations within the NJL model (see below). Therefore one 
can use the comparison of these two models to conclude about the importance of these 
corrections. 

After integrating eq. (§) by parts the contribution of valence quarks and antiquarks can 
be rewritten as 

Pciq{T, /i) = ^ ^ K(T, /x) + ng-(T, /i)] , (8) 

where and are the quark and antiquark occupation numbers, 

n.(T, /.) = ,^^^p[(^_^)/^] , -.-(T, /.) = -,(T, -;.) . (9) 
The baryon-chemical potential is determined by the net baryon density 

" ~3 9/i ^ 6n^ J ^^^^[^5*^^'^^ -nqiT,ix)]. (10) 



The net quark density is obviously n = Shb- The values for the a and vf fields and thereby 
the quark masses in eq. are obtained by minimizing fl with respect to a and n, 



on _ 

dQ 



v^)a-H + gps = , 



\2/ 2 , ^2 2\ I 



(11) 



(12) 



The scalar and pseudoscalar densities can be expressed as: 

d^p 1 



Ps = {qq) = gcriyq J j^^MT, p) + ng{T, p)] 



Pps = {qil^rq) = gnu^ 



d^p 1 



ng(T, p) +nq{T, fi)] 



(13) 



(14) 



The minima of Q defined by eqs. (|ll|), ( [12|) correspond to the stable or metastable states of 
matter in thermodynamical equilibrium where the pressure is P = —Qrnin- The a and pion 
masses are determined by the curvature of fl at the global minimum: 



Mi 



dirt 



2 • 



(15) 



Explicitly they are given by the expressions 
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2 '^Q 



27r2 
M2 , 



9^ 



dpp^ 



+ 



Vl + exp[(E + /i)/r] l + exp[(E-//)/T]^ 
1 1 



+ 



TE2 \^2(1 + cosh[(E + /i)/T]) 2(l + cosh[(£;-/x)/T]) 



(16) 



9^ 



27r^ 



dpp — 



1 + exp [(E + ^)/T] 



l + exp[(^-/i)/T] 
(17) 

2^2 



Here we have set the expectation value of the pion field to zero, vf = 0, thus M = g a . This 



version of the sigma model was used earlier in ref. for thermodynamical calculations at 
nonzero T and /x = 0, and at nonzero /i and T = 0. Some useful formulae for the case of 
small quark mass are given in the Appendix. 



III. THERMODYNAMICS OF THE NJL MODEL 



The NJL model has been widely used earlier for describing hadron properties and the 
chiral phase transition [T^|T^. The simplest version of the model including only scalar and 
pseudoscalar 4-fermion interaction terms is given by the Lagrangian 

G 



mo)q + — 



(qq) + {qi-i^rqY 



(18) 



where mo is the small current quark mass. At vanishing uiq this NJL Lagrangian is invariant 
under chiral SU{2)l x SU{2)fi transformations. The coupling constant G has dimension 
(energy)"^, which makes the theory non-renormalizable. Therefore, a 3-momentum cutoff 
A is introduced to regularize divergent integrals. It defines an upper energy limit for this 
effective theory. Free parameters of the model are fixed to reproduce correctly the vacuum 
values of the pion decay constant (93 MeV), pion mass (138 MeV), and the constituent 
quark mass (337 MeV). Below we use the following parameters [|18|: G = 5.496 GeV~^, 
mo = 5.5 MeV, and A = 631 MeV. With these parameters the chiral transition occurs at the 
temperature T ^ 190 MeV [|T^] (for /i = 0) which is significantly higher than in the sigma 
model. 

The partition function for the NJL model reads: 



Z = Tt exp — (H — fiAf 



VqVq exp 



(19) 



In the mean-field approximation the Lagrangian ( |18D is represented in a linearized form 



T|,[Ti: 



C = q{iYd^, - mo)q + G{qq){qq) 



G 



(20) 



such that the partition function becomes 



^As demonstrated in ref. |17| , the inclusion of the vector-axialvector terms may change signifi- 
cantly the parameters of the chiral phase transition, in particular, the position of the critical point. 
But this does not change the qualitative conclusions of the present paper. 



exp 



VG{qqf 
T 



detr 



{p,Y + ^^l'-M)/T 



(21) 



where the constituent quark mass is determined from the gap equation 

M = mo - G{qq) . 

The right-hand side of this equation involves the scalar density 

d^p 1 



(22) 



Ps 



(qq) = Mug J — -^MT, fi) + nq{T, /i) - 1] , 



(23) 



lp<A (27r)3 

where Ug and rig are the valence quark and antiquark occupation numbers defined in eq. (^. 
Here the last term in brackets gives the contribution from the Dirac sea (which corresponds 
to the vacuum part of the sigma model) and cannot be neglected. The rest comes from 
valence quarks and antiquarks similar to the sigma model (compare with eq. (|13[)). 
From eq. (|20|), the grand canonical potential for the NJL model can be written as: 



2G ''Vp<A(27r)3 



E + T\n 



1 + exp 



E + fi 



+ Tln 



1 + exp 



E- ^1 



T 



(24) 



The minimization of with respect to M gives the gap equation (|2^) . The expression (plf ) 
is formally identical with eq. (^ derived for the sigma model, but now the first term in 
curly brackets, coming from the Dirac sea, is treated explicitly after introducing the cut-off 
momentum A. One can calculate this vacuum contribution explicitly. 



p2 + M2 



UgA' 



v/rT^(H--J-^ln 



VTT^+ 1 



z 



(25) 



'Jp<A (27r)3 ■ - 87r2 

where z = M/A. Expanding this expression in powers of z one can find contributions to Q 
of order M^, M*^,.. . As mentioned above, in the sigma model the vacuum terms are partly 
absorbed in the coefficients of the effective potential U{(J, vf). However, the logarithmic term 
In cannot be removed in this way. Therefore, the NJL model has additional nonlinear 
terms in the vacuum energy which are responsible for the differences in the thermodynamic 
properties of the two models. 
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The sigma and pion masses are not as straightforward to obtain as in the hnear sigma 
model because in the NJL model they are not represented as dynamical fields. In this model 
mesons are described as collective qq excitations. Their masses can be obtained from the 
poles of the quark-antiquark scattering amplitude which can be computed, for instance, 
in Random Phase Approximation (RPA) In this way, one can derive the following 

equations for the sigma and pion masses, 

= ^ + (M,2 - AM')G liM, M.) , (26) 
The function I{x, y) is the quark-antiquark propagator defined as: 




where E = y/x'^ + p'^, and the occupation numbers defined in (||). In this 

integral V means principal value. 



IV. NUMERICAL RESULTS 

In the case of the linear sigma model everything is determined when the gap equation (pUD 
is solved in the (T, fi) plane, whereas in the NJL model the meson masses have to be solved 
for as well. Below we present results of our numerical calculations (see also Appendix). 



A. Phase diagrams 

We start this section with presenting in Fig. |1] the resulting phase diagrams in the (T, fi) 
plane calculated for the two models. The middle line corresponds to the states where the 
two phases co-exist in the first order phase transition. Along this line the thermodynamical 
potential Q has two minima of equal depth separated by a potential barrier which height 
grows towards lower temperatures. At the critical point C the barrier disappears and the 
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transition is of second order. The other hnes in Fig. |l| are spinodal hnes which constrain 



this instabihty can be obtained from dynamical simulations @J^. 

It is instructive to plot the thermodynamic potential as a function of the order parameter 
for various values of T at = 0, and for various values of yU at T = 0. The first case is 
shown in Fig. ^, where the left panel is for the sigma model and the right panel for the 
NJL model. One clearly sees the smooth crossover of the symmetry breaking pattern in 
both cases. Note that the effective bag constant (the energy difference between the global 
minimum and the local maximum of the potential in vacuum) is about 100 MeV/fm^ in the 
NJL model, whereas in the sigma model it is significantly smaller, ~ 60 MeV/fm^. To a large 
extent this difference is responsible for the difference in the temperatures corresponding to 
the crossover transition: about 140 MeV in the sigma model and about 180-190 MeV in the 
NJL model. 

In Fig. ^the same plot is shown for T = and a nonzero fi. Here, one clearly observes the 
pattern characteristic for a first order phase transition: two minima corresponding to phases 
of restored and broken symmetry separated by a potential barrier. The barrier height is 
larger in the sigma model than in the NJL model, thus, indicating a weaker first order phase 
transition in the NJL model. It now follows that somewhere in between these two extremes, 
for some fic and Tc, there exists a second order phase transition (the critical point). Indeed, 
this point is found and shown in Fig. 0. The corresponding values are (Tc, fic) — (99, 207) 
MeV in the sigma model, and (Tc, fic) — (46, 332) MeV in the NJL model. The behavior of 
the thermodynamic potential at fi = fic and various T is shown in Fig. ^ One can see that 
the potential has only one minimum which is flattest at the critical point. 



Now let us consider the model predictions for the effective masses. The constituent quark 



the regions of spinodal instability where 




T 



< 0. Information about the timescales of 



B. Effective Masses 



mass is shown in Fig. ^ as function of T and /x. These plots, of course, show the same phase 
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structure as discussed above. At /i = in both models the quark mass falls smoothly from 
the respective vacuum value and approaches zero as T goes to infinity. One could define a 
crossover temperature as corresponding to a steepest descent region in the variation of M. 
This again gives a temperature of about 140 — 150 MeV for the sigma model and about 
180 — 190 MeV for the NJL model. At T = and nonzero /i the constituent quark mass 
shows a discontinuous behavior reflecting a first order chiral transition. 

The sigma and pion masses for various T and are shown in Figs. In both models 
the sigma mass first decreases smoothly, then rebounds and grows again at high T. The pion 
mass does not change much at temperatures below Tc but then increases rapidly, approaching 
the sigma mass and signaling the restoration of chiral symmetry. As T goes to infinity the 
masses grow linearly with T. The fi = fic case is especially interesting in the sigma model. 
Since the sigma field is the order parameter of the chiral phase transition, its mass must 
vanish at the critical point for a second order phase transition. This means that Q has zero 
curvature at this point. It is, however, not clear what the sigma mass should be at the 
critical point in the NJL model where the quark condensate (qq) is the order parameter. 
Fig. ^ indeed shows that exactly at the critical point the sigma mass is zero in the sigma 
model. This is not the case in the NJL model, at least within the RPA used here. 

In Fig. 1^ the masses are plotted as function of /i for T = and T = Tc. For T = 
one clearly sees discontinuities in the behavior of the masses characteristic for the first order 
phase transition. 

An interesting point is that, in the linear sigma model, there is no stable phase with 
heavy quarks for T = 0, i.e., the quark mass assumes its vacuum value all the way up to 
the chiral transition, and then drops to a small value in the phase where chiral symmetry 
is restored (see Fig. |^. This behavior is related to the appearance of a bound state at 
zero pressure. Within the linear sigma model this "abnormal" bound state was found by 
Lee and Wick a long time ago |^ . Recently, it was shown in ref. [|17| that a similar bound 
state appears also in the NJL model. This behavior, however, depends on the value of the 
coupling constant g or G. For our choice of g and G, this state exists in the linear sigma 
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model, but not in the NJL model, where there is a stable phase of heavy quarks at T = 0, cf. 
Figs. ^ and |^. In general, if the coupling constant is sufficiently large, then the attractive 
force between the constituent quarks becomes large enough to counterbalance the Fermi 
pressure, thus giving rise to a bound state. To demonstrate this we have varied the coupling 
constant g for the sigma model within reasonable limits. 

The results for the quark mass are shown in Fig. ^ It is seen that, indeed, one can 
change the smooth crossover for /i = into a first order transition by increasing the coupling 
constant (Fig. || left panel) and change the first order transition in the case of T = into 
a smooth crossover. In this way a heavy quark phase comes into existence as the coupling 
constant is decreased and the bound state disappears (Fig.|| right panel). An analogous 
investigation for the NJL model leads to similar results. 



C. Adiabats 

Regarding hydrodynamical simulations the entropy per baryon is an interesting quantity. 
One can easily calculate it using standard thermodjTiamic relations, 

A Tn ' ^ ^ 

where e, p, n are respectively the energy density, pressure and net density of the quarks and 
antiquarks {n = Sub). By studying this quantity, one can check if there is a tendency 
towards convergence of the adiabats towards the critical point as was claimed in ref. If 
this was the case it would be easy to actually hit or go close to this point in a hydrodynamical 
evolution. Fig. ^ shows the contours of S/A in the (T, /i) plane calculated in the sigma model 
(left) and in the NJL model (right). We actually observe a trend which is quite opposite to 
this expectation. It turns out that the adiabats turn away from the critical point when they 
hit the first order transition line and bend towards the critical point only when they come 
from the smooth crossover region. This is explained as follows. First, note that all adiabats 
terminate at zero temperature and /x = Myac, i-e. the (T, /z) combination corresponding to 
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the vacuum. The reason is that as T — 0, also 5* ^ (by the third law of thermodynamics), 
therefore, for fixed S/A we have to require that n — > 0, which is fulfilled when fi = M^ac- 
For our choice of parameters, in the sigma model the point {T,fi) = (0,M^ac) is also the 
endpoint of the phase transition curve at T = 0, since the phase transition connects the 
vacuum directly with the phase of restored chiral symmetry, cf. Figs. ^ and ^. For the 
NJL model, the endpoint of the phase transition curve is not identical with (0, M^ac) but is 
rather close to it. Therefore, also the adiabats which hit the phase transition curve have to 
bend away from the critical point and approach the endpoint of the phase transition line at 
T = 0, i.e., T decreases and /i increases. 

This behavior is quite opposite to the case underlying the claim in ref. 0, where the 
hadronization of a large number of quark and gluon degrees of freedom into relatively few 
pion degrees of freedom leads to the release of latent heat and consequently to a reheating 
(increase of T) through the phase transition. Remember, however, that in our case there 
is actually no change in the number of degrees of freedom in the two phases, only in their 
respective masses. Consequently, there is no "focusing" effect in the linear sigma and NJL 
models. 

V. CONCLUSIONS 

We investigated the thermodynamics of the chiral phase transition within the linear 
sigma model coupled to quarks and the Nambu-Jona-Lasinio model. These models have 
similar vacuum properties but treat the contribution of the Dirac Sea differently. In the 
sigma model this contribution is "renormalized out" while in the NJL model it is included 
explicitly up to a momentum cutoff A. By comparing thermodynamic properties of these 
two models one can check the importance of these vacuum terms. In both models, we 
found for small bare quark masses a smooth crossover for nonzero temperature and zero 
chemical potential and a first order transition for zero temperature and nonzero chemical 
potential. The first order phase transition line in the (T, /z) plane ended in the expected 
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critical point. It has been found that the a mass is zero at the critical point in the sigma 
model whereas in the NJL model it always remains nonzero. The phase transition in the 
sigma model turned out to be of the liquid-gas type. This, however, depends on the coupling 
constant g between the quarks and the chiral fields. Prom the comparison we conclude that 
the phase transition pattern is generally weaker in the NJL model than the sigma model. 
Certainly, it will be interesting to use both models in hydrodynamical simulations in order 
to confirm or disconfirm possible observable signatures of the phase transition discussed in 
the introduction. In particular, the sigma model which contains dynamical a and pion fields, 
would be suitable to study the long wavelength enhancement of the a field at the critical 
point. Such simulations are in progress. 
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In the chirally symmetric phase the constituent quark (antiquark) mass M is small, it 
goes to zero in the chiral limit. Therefore, it is instructive to evaluate the thermodynamic 
potential 0(T, /x, M) for small M. In this limit Vtqq can be represented as a power series in 
M. Below we give explicit expressions for the sigma model. Taking into account that is 
an even function of M, one can write 
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APPENDIX A: 




(Al) 
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Here the first term, flo(T,fi) = Qgg{T, ^,0), can be easily calculated for arbitrary T and /i. 
The well-known result is 



27r2 



(A2) 



180 6 ^ 12' 

The quark number and entropy densities for massless fermions are obtained by differentiating 
flo{T,fj,) with respect to fj. and T, respectively, 

n = -p^(vr^TV + /i='), (A3) 



Ug ( 7tt' 



4 



-T^ + n^Tfx^ . (A4) 



Gtt^ \ 15 

The second term in eq. ( |Al| ) differs only by a factor of M from the scalar density defined 
in eq. (pTSI). A straightforward calculation gives 



This can be used to estimate the pion and sigma masses at large T and/or /i. Expressing 
in terms of mean it and a fields, eq. ([T5|) , and using the definition of effective masses 
from eq. (^, one arrives at the following asymptotic (M — > 0) expression for the pion and 
sigma masses 

Mi = M! = g\ + ^) . (A6) 

It shows that deep in the chiral symmetric phase the pion and sigma masses are degenerate 
and large. At high temperatures (T ^ /i), M^^ = = gT, where ~ 3 in our calculations. 
Therefore, the contribution of pion and sigma excitations to the thermodynamical potential 
is negligible. 

In case of the NJL model the above expressions are slightly modified due to the finite 
cut-off A in the momentum integration. 
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FIG. 1. The phase diagrams for the sigma model (left) and the NJL model (right) in the (/x,r) 
plane. The middle curve is the critical line and the outer lines are the lower and upper spinodal 
lines. C is the critical point. 
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FIG. 2. The thermodynamical potentials Q for the sigma model (left) and the NJL model 
(right). For both models jj. = 0. The levels correspond to (starting from the top): 

T = [0, 100, 135, 155, 175, 190] MeV for the sigma model and T = [0, 100, 140, 170, 200, 230] McV 
for the NJL model. 
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FIG. 3. The thermodynamical potentials Q for the sigma model (left) and the NJL model 
(right). For both models T = 0. The levels correspond to (starting from the top): 
H = [0, 225, 279, 306, 322, 345, 375] MeV in the sigma model, and /x = [0, 288, 343, 348, 35, 378, 396] 
MeV for the NJL model. 
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FIG. 4. The thermodynamical potentials J7 for the sigma model (left) and the NJL model (right). 
For the sigma model case jj, is fixed to 207 MeV and the levels correspond to (starting from the 
top): T = [0, 50, 75, 100, 125, 150] MeV. For the NJL model case n is fixed to 332 MeV and the 
levels correspond to (starting from the top): T = [0, 28, 46, 70, 105, 133] MeV. 




FIG. 5. The constituent quark (antiquark) mass in the sigma model (left) and the NJL model 
(right) as function of n and T. 
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FIG. 6. The sigma mass (solid line) and pion mass (dashed line) in the sigma model (left) and 

NJL model (right) as functions of temperature for /x = (right pair) and for /x = /Xc (left pair). 




FIG. 7. The sigma mass (solid line) and pion mass (dashed line) in the sigma model (left) and 



NJL model (right) as functions of chemical potential for T = (right pair) and for T = Tc (left 
pair). 
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FIG. 8. The constituent quark (antiquark) mass as function of temperature for zero chemical 
potential (left) and as function of chemical potential (right) for zero temperature in the sigma 
model. The solid line represents the mass for g = 4.5, the dashed-dotted line for g = 3.3 and the 
dashed line for ^ = 2.8. 
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FIG. 9. The entropy per baryon number, S/A for the sigma model (left) and the NJL model 

(right). In the sigma model the curves correspond to (from left) S'/A=[28,21,17,13,ll,9,6,2]. In 

the NJL model the curves correspond to (from left) >S'/^=[22,16,11,8,6,5,3,1]. 
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